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1 INTRODUCTION. 



Picture three points at the vertices of an equilateral triangle in two dimen- 
sions, or four points at the vertices of a regular tetrahedron in three di- 
mensions. Thought of as scatterings of data they wouldn't seem to reveal 
strong linear associations between the coordinates. There are no clear axes 
of elongation in the scatterplots, which would suggest that change in some 
variable is predictable as a function of the others. In general, such associa- 
tions are usually indicated by the covariance matrix S u of the set of points 
u = {xi, . . . , x n } in R p , which is given by 



The off-diagonal entries of S u , the pairwise covariances, tell us something 
about dependencies. If the coordinate variables are independent these en- 
tries are zero. Though the converse is false, a diagonal covariance matrix 
roughly says that the coordinates are not mutually linearly predictable from 
each other. Indeed, for our equilateral triangle in 1R 2 , tetrahedron in R 3 , 
and the generalized configurations in higher dimensions having equal inter- 
point distances, the covariance matrix turns out to be diagonal. In fact, it's 
a scalar multiple of the identity. Furthermore, the converse is also true: any 
configuration of n = p + 1 points in p dimensions whose covariance matrix 




where 
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is a positive multiple of the identity are equidistant from each other. We 
formalize this result in the following theorem: 

Theorem. Let u = {x 1; . . . , x n } be a set ofn points in M p , with n = p + 1 > 
2, and let a 2 be an arbitrary positive number. Then the interpoint distances 
of u satisfy ||xj — Xj|| 2 = 2o~ 2 5ij if and only if nS u = a 2 l p . 

In other words, p+1 points in p dimensions lie at the vertices of a regular 
simplex if and only if their covariance matrix is a multiple of the identity. 
A proof of this statistical characterization of regular simplices is given in 
section |2j after some preliminaries. 

2 STATISTICAL CHARACTERIZATION OF 
REGULAR SIMPLICES. 

The reader is assumed to be familiar with the basic elements of linear al- 
gebra in R p (linear subspaces, span, linear dependence and independence, 
basis and dimension), as treated, for example, in the text of Seber pp. For a 
finite subset u of R p let V u = span{x — x u : x e u}. 

Lemma. With n > 1 let u be any collection of n points in M p with common 
squared interpoint distance 2a 2 > 0. Then dim(V u ) = n — 1, and with 
r 2 n = a 2 {n — l)/n and s 2 n = a 2 jinin — 1)), the following are true for each 
x in u; 

I |x X u | | ^gn, | |x u X u _| x | 1 1 Su^n, X X u _j x j _L V u _j x }.. 

Proof. We argue by induction. The three claims are easily verified if n — 2. 
When n > 2, for every x in u the points u — {x} are equidistant from x, and 
by the induction hypotheses also equidistant from their average x u _{ x }, albeit 
at a smaller distance. Hence, the points of u — {x} lie on the intersection of 
two spheres with distinct centers, x and x u _{ x }, which implies that V u _{ x } is 
perpendicular to the direction vector of the line 

L u , x ( a ) = x u-{x} + at (x - x u _ {x} ) (a E R) 
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passing through these centers and that the points of u — {x} are equidistant 
from each point of L u x . In particular, all points of u — {x} are equidistant 
from L U)X (l/n) = x u , hence so are all points of u. Because x — x u _{ x } _L 
V u -{ x } but x — x u _{ x } G V u , dim (V u ) = dim (V u _{ x j) + 1. By orthogonality 
Il x u — x u-{x}|| 2 = r ln ~ r ln-i an d does not depend on x. Using the fact 
that x u _{ x },x u , and x all lie on L ux in tandem with orthogonality gives 
2cr 2 = r 2 n _ 1 + (s a>n + r^ n ) 2 ] solving these two equations for r^ n and s a>n 
finishes the induction. ■ 

Proof of the theorem. Let X = (xi, . . . , x n ), an element of R pxn . Since 
St(u) — S u for any translation T, we can assume without loss of generality 
that the members of u have already been centered by subtraction of their 
mean, so x u = and in general letting B v := |v|S v we have 

B u = ^xx' = XX'. (1) 

x£u 

Assuming that the points are equidistant, we infer from (CQ) and the lemma 
that 

b u x = yy' x + xx ' x = ( r ln -° 2 ) y + r l, n * 

yeu-{x} yeu-{x} 
I 2 2 \ , 2 2 

for each x in u. Hence B u x = a 2 I p x on V u . Since dim(V u ) = p by the 
lemma, B u = cr 2 I p . 

For the converse, assume that B u = cr 2 I p . Note that the matrix 

A = a" 2 X'X G R nxn 

is symmetric, A' = A, and idempotent, A 2 = cr" 4 X'XX'X = o- _4 X'B u X = 
A. Hence A is an orthogonal projection, and therefore has rank equal to its 
trace, 

rank(A) = tr(A) = (j- 2 tr(XX) = a- 2 tr(XX) = a~ 2 tr(B u ) = p, 

using the cyclic invariance of the trace. With l n G R n the vector with all 
components equal to 1, Al„ = by virtue of x u = 0. By the rank plus 
nullity theorem the null space of A has dimension one, and must therefore 
equal span(l n ), the span of l n . Hence A = I„ — - l n l' n , as this is the unique 
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orthogonal projection of rank p with null space span(l„). As the entries of A 
are a~ 2 times the inner products of the vectors in u, the squared interpoint 
distances between Xj 7^ Xj equals 

||x l -x j || 2 = 2(x^x i -x i x,) = 2(T 2 N + =2a 2 . 



We remark that once the matrix A is determined to have constant off- 
diagonal entries, the proof may also be completed by induction in the fol- 
lowing more geometric way: Assume that n > 2, the base case being trivial. 
Any p points in IR P lie in a hyperplane of dimension p—1, and for x G u let 
H denote the hyperplane which contains u — {x}, the space V u _{ x } trans- 
lated by x u _{ x }. The inner products x'y for all y G u — {x}, being the 
off-diagonal elements of A, are equal, and therefore, x'y = x'x u _{ x }, so 
x _L y — x u _{ x }. Hence x _L V u _{ x }, and since x u _{ x } = — x/p, we conclude 
x — x u _{ x } _L V u _{ x } = H — x u _{ x }. 

Now let T be the translation Ty = y — x u _{ x }, and, with {ej}i<j< p the 
standard basis, O the rotation that maps Tx to (3e p where (3 = | |x — x u _{ x }| |. 
That is, V(x) = (3e p for V = OT, and 

B V (u-{x}) + P 2 e p e' p = O (B T(u _ {x)} + Tx(Tx)') O' (2) 
= O (B u {x} + Tx(Tx)') O' = OB T( u)0' = OB u O' = a 2 OI p O' = a%. 

Since V u -{ x } _L x — x u _{ x }, V(H) C BP' 1 x {0}, and we can consider the 
points V(u — {x}) as lying in H p_1 . By (j2J), the (p — 1) x (p — 1) submatrix 
[Bv( u -{x})]i<ij<p-i equals <j 2 I p _i, so applying the induction hypotheses to 
V(u— {x}) we conclude that the interpoint distances of u— {x}, unchanged 
by V, are all 2<r 2 . The induction is completed by noting that this is true for 
each x in u. 
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